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By way of a complete integration of the Friedmann equations, in terms of observables, it is shown 
that for the cosmological constant A > there exist non-flat FLRW models for which the total 
density parameter Q remains ~ 1 throughout the entire history of the universe. Further, it is 
shown that in a precise quantitative sense these models are not finely tuned. When observations 
are brought to bear on the theory, and in particular the WMAP observations, they confirm that we 
live in just such a universe. The conclusion holds when the classical notion of A is extended to dark 
energy. 



The flatness problem is often considered to be the most 
impressive issue in standard cosmology that is addressed 
by the inflation paradigm p). Let us start by summariz- 
ing the flatness problem for zero cosmological constant. 
With A = the density parameter of a FLRW model is 
given by 



n = 



8irGp 



(1) 



and for a single fluid the state space is summarized for 
standard models in FIG. ^ 3- 
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FIG. 1: The state space for a single fluid FLRW model with 
p = (7— l)p and 7 > 2/3. M is the Milne universe (Minkowski 
space). Bb is the spatially flat FLRW model and the Big 
Bang. Turning points (H — 0) are sent to infinity by the def- 
inition Q. The diagram reflects the evolution for expanding 
universes. 

The essential point is that except for the spatially flat 



case 



fi = fi(i). 



Observations show that 



fi ~ 1. 



(2) 



(3) 



The flatness problem involves the explanation of given 
(J2J. The problem can be viewed in two ways. First, one 
can take the view that there is a tuning problem in the 
sense that at early times fi must be finely tuned to 1 Q . 
However, this argument is not entirely convincing since 
all standard models necessarily start with fi exactly 1. 
More convincing is the view that except for the spatially 
flat case the probability that fi ~ 1 is strongly dependent 
on the time of observation [3 and so there is an epoch 
problem: why should © hold? 



In this letter we point out that if A > then there 
exist standard models for which fi ~ 1 throughout their 
entire evolution even though they are not spatially flat. 
Moreover, in a precise quantitative sense, we show that 
these models are not finely tuned. When current obser- 
vations are brought to bear on the theory, they confirm 
that we live in just such a universe [||. 



To include A with dust define, in the usual way, 
c 2 A „ c 2 k 



fi ' 



ill 



3H 2 ' " K ~ H 2 R 2 
so that the Friedmann equations reduce to 

n M + n A + n k = i 



(4) 



(5) 



where we have written fi = Q,m f° r convenience. To 
complement FIG.^ which is the state space for ft\ = 0, 
the state space for J7m = is shown in FIG. [3 
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FIG. 2: The state space for f!jn = (here rotated left n/2 to 
save space). Again M is the Milne universe (Minkowski space 
0)- dS is the spatially flat representation of de Sitter space 
and acts as an attractor for the two other representations of 
de Sitter space. 

Evolution in the full f2^_ — VLm plane was first considered 
by Stabell and Refsdal |8| for the case of dust (7 = 1). 
Writing Z = 1 + z (1 + z the standard redshift), they 
populated the phase plane from relations equivalent to 



fi 



M 



n Ma z 3 + (1 - fi Mo - fiAjz 2 + fiA 



(6) 



which is a restatement of the conservation law for dust 
(the constancy of VImH 2 R z ), and 



fi A 



fi/ 



fi Mo Z 3 + (1 - fi 



Mo 



fi A )z 2 + n r 



(7) 



which is a restatement of the constancy of A. They dis- 
tinguished trajectories via (half) the associated absolute 
horizon, J °° The same technique has been applied 
more recently and in more general situations 0. 

Equivalcntly, from a dynamical systems point of view 
jlCj . we can, again using dust as an example, consider 
the system of differential equations 



n A = (o M - 2n A + 2)n A , 



and 



M 



(n 



20/ 



1)0 



M ■ 



(8) 



(9) 



where = d/drj and r\ = ln(X/Z). We now recognize the 
critical points in FIG.n an d FIG. [21 Bb is a repulsor, dS 
an attractor and M a saddle point for the full Oa — CIm 
plane. 

The approach used here is somewhat different. Al- 
though what now follows can be generalized ^1 > we con ~ 
tinue with dust as it presents an uncluttered relevant ex- 
ample. We observe a constant of the motion (a) for the 
system (JSJ-I©. The constant is given by |l2| 



( 



207 } ±a{ 



30a" 



) 3 = o 



(10) 



for k = ±1 respectively. Since Oa = 1 — O^/ for k = 0, 
the system can be considered solved and the Friedmann 
equations, in terms of observables, in effect integrated. 
A phase portrait with trajectories distinguished by a is 
shown in FIG.EJ 

The physical meaning of a is not hard to find. In terms 
of the bare Friedmann equation 



R 2 + k = 



C AR 2 



R 



(11) 



with (k, A) considered given, each value of the constant 
C > determines a unique expanding universe (k =/= 0). 
For each C, if k = 1, there is a special value of A that 
gives a static solution (the Einstein static universe or 
asymptotic thereto), 



A E = 



9C 2 ' 



In terms of A^ we have 



A 

AT 



(12) 



(13) 



and so a is a measure of A relative to the "Einstein" 
value The state spaces shown in FlC^and FIG. [21 
correspond to a = 0: A = in the first case and C = 
in the second. 

As FIG. makes clear, if A > and k ^ then 

n = o M + n A ~i (14) 




FIG. 3: The loci a — constant and Ht = constant [l3^ | 
in the — Q.m plane. The loci shown are (bottom up) 
Ht = 0.6, 2/3, 3/4, 0.826, 1, 1.5, oo (= a = 1) 
(thick) and \a\ = 0.01, 0.05, 0.1, 1/4, 1/2, 1 (thick), 
2, 3, 4, 6, 8, 12, 20, 40, 100, 500. Note that for A > 
the delimiter Qa = 1 — Am (k — 0) is approached on either 
side (k — ±1) as a — > oo. 



throughout the entire evolution of the associated universe 
if a G (a, oo) where a is a matter of choice but of the order 
>^ 500. In any event, it is clear that a need not be finely 
tuned to produce l|14|) 15]. 

Fortunate we are to live in an era of unprecedented 
advances in observational cosmology. We now bring 
some of theses observations to bear on the foregoing dis- 
cussion. The Wilkinson Microwave Anisotropy Probe 
( WMAP) has enabled accurate testing of cosmological 
models based on anisotropics of the background radia- 
tion [16|]. Independently, the recent Hubble Space Tele- 
scope (HST) type la supernova observations [I3 not only 
confirm earlier reports that we live in an accelerating uni- 
verse 18] 19], but also explicitly sample the transition 
from deceleration to acceleration. A comparison of these 
results in a partial phase plane is shown in FIG. 0] 20]. 
Clearly the ( WMAP) results, and to a somewhat lesser 
extent the (HST) results [2l|, support the view that we 
live in "large" a universe for which l|14l) has held through- 
out its entire evolution. 

We have argued that there are an infinity of standard 
non-flat FLRW dust models for which O ~ 1 through- 
out their entire evolution as long as A > 0. Further, we 
have shown that the WMAP observations confirm that 
we live in just such a universe. The idea can be gen- 
eralized in a straightforward way to more sophisticated 
multi-component models wherein CI ~ 1 (fi = for 
all species i but not flk) from Bb to dS within a phase 
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FIG. 4: WMAP and HST results superimposed on 
part of the phase portrait of FIG. [3] In the former case the 
"WMAP only" results are used and in the latter case earlier 
confidence levels have been removed for clarity. The values of 
Hoto shown are (top down) 1.5, 1, 0.826, 3/4 and 2/3. Above 
the delimiter £Ia = 1 — !)jj , k = 1 and the values of a shown 
are (bottom up) 500, 100, 40, 20, 12, 8, 6 and 5.4. Below the 
delimiter k = — 1 and the values of a shown are (top down) 
500, 100, 40, 20, 12, 8, 6 and 5.4. 

portrait hyper-tube about k = (see 0] for a two- 
component model). Note, however, that the dust model 



examined here is an excellent approximation now and 
the model relevant to a comparison with current obser- 
vations. With Q ~ 1 there is no tuning or epoch prob- 
lem and so no flatness problem in the traditional sense. 
However, our presence in this hyper-tube is presumably 
favored and an explanation of this probability about the 
k = trajectory in a sense presents a refinement of the 
classical flatness problem. 

Since the analysis here has been based on the classical 
notion of A, it is appropriate to conclude with a query as 
to whether or not this analysis is stable under perturba- 
tions in the definition of A itself (that is "dark energy" 
as opposed to the classical notion of A). Since A can 
be introduced by way of a component p w = wp w where 
w = —1, consider w unspecified (but < — 1/3). It can 
then be shown that the analysis given here generalizes 
in a straightforward way |22| and that for perturbations 
about w — — 1 the phase portrait given in FIG. [31 is sta- 
ble, but the state space for Q,m = shown in FIG. |3 
is not 23]. This latter point in no way alters the fact 
that there remain an infinity of non-flat FLRW models 
for which Q ~ 1 throughout their entire evolution as long 
as W{= 8irp w R 3( - 1+w ^) > 0. 
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Introducing the (non-interacting) species p w = wp m (w < 
— 1/3) with the dust and with no concern here for the 



classical energy conditions, we have the constant W = 
S-Kp w R 3{l+w) and so define fi w = c 2 W /ZH 2 R 3{1+w) to 
replace Oa- It follows that equations JSJ, (7J, JHJ and © 
are replaced by 

Mm z 3 + (i - n Mo - n Wo )z 2 + nwcZ^+v) ' 



Q. Wn z' M1+w) 



[23] 



Mm Z 3 + (1 - n Mo - ttwJZ 2 + Q. Wo Z^+™) ' 



tt w = (Sl M + (1 + 3w)Q w - (1 + 3w))Q w , 



and 



(« M + (i + 3w)n w - i)a, 



respectively. More importantly, llHH generalizes to 



g, Q.M |-i-3«i _^ a , 3m . 
Qw w 1 + 3w 



i+3w f \3w _ q 



for k — ±1 respectively. Now a — W/We where 
W E = --(- 3wC 



w 1 + 3w 

The stability of the phase portrait under perturbations 
on w now follows directly. 

The instability of the state space for Qm = 0, and there- 
fore the attractor for the full phase plane, is not difficult 
to see. With Qm = k = the Friedmann equation is 



W 







where e = —3(1 + w) and here we take W > 0. This 
equation is unstable wrt e in the sense that the qualitative 
behaviour of the solutions changes with e. Whereas we 
obtain the familiar exponential solution of de Sitter space 
for e = 0, for t/Owe obtain 



where c is a constant. The associated spacetimes are sin- 
gular at t = c ( for example, the Ricci scalar diverges like 



(*- 



For expanding solutions, with e>0(ui<— 1) 



t G ( — co, c). The universe is asymptotically past flat but 
ends in a "Big Rip" (bR) at t = c. (See R. C. Caldwell, 
M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 
91, 071301 arXiv:astro-ph/0302506 .) We replace dS 
by bR in the associated phase plane and observe that all 
standard universes end at bR. The inclusion of a classi- 
cal A does not change this. For expanding solutions with 
e < (w > —1) t 6 (c, co). The universe starts at a 
big bang but is future asymptotically flat. In this case 
the inclusion of a classical A reverts the attractor to dS. 
It can be shown that these conclusions are unchanged 
if w = w(t). Whereas it is clear that observations can 
never prove that k = (but they can prove k ^ 0), it is 
also clear that observations also can never determine the 
asymptotic future state of the universe. 



